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Abstract 

We investigate an original domain- wall model in 4+1 dimensions numerically 
in the presence of U(l) dynamical gauge field only in an extra dimension, 
corresponding to a weak coupling limit of 4-dimensional physical gauge cou- 
pling. Using a quenched approximation we carry out numerical simulation 
for this model at f3 s {= l/<7s) = 0.29 ("symmetric" phase) and 0.5 ("broken" 
phase), where g s is the gauge coupling constant of the extra dimension. In 
the broken phase, we found that there exists a critical value of a domain-wall 
mass mjj which separates a region with a fermionic zero mode on the domain 
wall from the one without it in the same case of (2+l)-dimensional model. 
On the other hand, in the symmetric phase, our numerical data suggest that 
the chiral zero modes disappear in the infinite limit of 4-dimensional volume. 
From these results it seems difficult to construct the U(l) lattice chiral gauge 

theory via an original domain- wall formulation. 
11.15Ha, 11.30Rd, 11.90. +t 
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I. INTRODUCTION 



Although perturbative aspects of electro-weak interaction are well described by the stan- 
dard model, its non-perturbative phenomenon such as baryon number asymmetry have to be 
investigated beyond the perturbation theory. At present, the most powerful non-perturbative 
technique is the lattice field theory. However it is non-trivial to define the standard model 
on a lattice since it is a type of chiral gauge theories, construction of which is one of the 
long-standing problem of lattice field theory: Because of the fermion doubling problems, 
a naively D- dimensional discretized lattice fermion field yields 2 D fermion particles, half 
of one chirality and half of the other, so that the theory becomes non-chiral 0. Several 
lattice approaches have been proposed, but so far none of them have been proven to work 
successfully @. 

D. Kaplan has proposed a new construction of lattice chiral gauge theories via domain- 
wall models |J . Starting from a vector-like gauge theory in 2k + 1 dimensions with a fermion 
mass term being a shape of a domain-wall in the (extra) (2k + l)th dimension, he showed in 
the weak gauge coupling limit that a massless chiral state arises as a zero mode bound to the 
2A;-dimensional domain-wall while all the doublers have large masses of the lattice cut-off 
scale. It has been also shown that the model works well for smooth background gauge fields 



Two simplified variants of the original Kaplan's domain-wall model have been proposed: 
an "overlap formula" |§|| and a "waveguide model" . Gauge fields appeared in these 



variants are 2 k- dimensional and are independent of the extra (2k + l)th coordinate, while 
those in the original model are (2k + l)-dimensional and depend on the extra (2k + l)th 
coordinate. These variants work successfully for smooth background gauge fields [p^ - [T8| , |2T| , 
as the original one does. Non-perturbative investigations for these variants seems easier than 
that for the original model due to the simpler structure of gauge fields. 

However it has been reported |10|Jll| that the waveguide model in the weak gauge cou- 
pling limit cannot produce chiral zero modes needed to construct chiral gauge theories. In 



2 



this limit, if gauge invariance were maintained, pure gauge field configurations equivalent to 
the unity by gauge transformation would dominate and gauge fields would become smooth. 
In the set-up of the waveguide model, however, 2fc-dimensional gauge fields are non-zero 
only in the layers near domain wall (waveguide) , so that the gauge invariance is broken in 
the edge of the waveguide. Therefore, even in the weak gauge coupling limit, gauge fields are 
no more smooth and become very "rough", due to the gauge degrees of freedom appeared 
to be dynamical in this edge. As a result of the rough gauge dynamics, a new chiral zero 
mode with the opposite chirality to the original zero mode on the domain wall appears in the 
edge, so that the fermionic spectrum inside the waveguide becomes vector-like. It has been 



claimed [10,11] that this "rough gauge problem" also exists in the overlap formula since the 



gauge invariance is broken by the boundary condition at the infinity of the extra dimension 



17,18]. Furthermore an equivalence between the waveguide model and the overlap formula 



has been pointed out for the special case [19 . Although the claimed equivalence has been 



challenged in Refs. [^,^1|], it is still crucial for the success of the overlap formula to solve 
the "rough gauge problem" and to show the existence of a chiral zero mode in the weak 
gauge coupling limit. 

In the original model there are two inverse gauge couplings (3 = 1/g 2 and (3 S = l/g^, 
where g is the coupling constant in (physical) 2k dimensions and g s is the one in the (extra) 
(2k + l)th dimension. Very little are known about this model except (3 S = case |T0| , |22| j2! 



where the spectrum seems vector-like and the case of (2+l)-dimensional U(l) model p4 
Since perturbation theory for the physical gauge coupling g is expected to hold, the fermion 
spectrum of the model can be determined in the limit that g — > 0. In this weak coupling 
limit, all gauge fields in the physical dimensions can be gauged away, while the gauge field 
in the extra dimension is still dynamical and its dynamics is controlled by (3 S . Instead of the 
gauge degrees of freedom in the edge of the waveguide, (2k + l)th component of gauge fields 
represent roughness of 2/c-dimensional gauge fields. An important question is whether the 
chiral zero mode on the domain wall survives in the presence of this rough dynamics. The 
dynamics of the gauge field in this limit is equivalent to 2/c-dimensional scalar model with 



2L S independent copies where 2L S is the number of sites in the extra dimension. In general at 
large j3 s such a system is in a "broken" phase where some global symmetry is spontaneously 
broken, while at small f3 s the system is in a "symmetric" phase. Therefore there exists a 
critical point and it is likely that the phase transition at (3 S = /3° is continuous (second 
or higher order). The "gauge field" becomes rougher and rougher at smaller f3 s . Indeed we 



know that the zero mode disappears at j3 s = p2| , while the zero mode exists at f3 s = oo 
( free case ). So far we do not know the fate of the chiral zero mode in the intermediate 
range of the coupling f3 s . There are the following three possibilities: (a) The chiral zero 
mode always exists except /3 a — 0. In this case we may likely construct a lattice chiral gauge 
theory in both broken (J3 a > /3f) and symmetric (J3 S < f3f), and the continuum limits may be 
taken at j3 s = (3 C S . This is the best case for the domain-wall model, (b) The chiral zero mode 
exists only in the broken phase (f3 s > fif) . In this case the domain- wall method can describe 
a lattice chiral gauge theory in the broken phase at finite cut-off. However it is likely that 
the continuum limit taken at (3 S = (3 C S from above leads to a vector gauge theory, (c) No 
chiral zero mode survives except (3 S = oo. The original model can not describe lattice chiral 
gauge theories at all. It is very important to determine which possibility is indeed realized 
in the domain-wall model. 

Instead of (4+1) -dimensional models, we have recently investigated a (2+l)-dimensional 
U(l) model ||24|| . Using a quenched approximation we have carried out a numerical simulation 
to see whether chiral zero modes exist or not in this model. In the weak coupling limit of 
the physical gauge coupling the (2+l)-dimensional U(l) gauge system is reduced to the 
2-dimensional U(l) spin system []. Strictly speaking, there is no order parameter in the 
2-dimensional U(l) spin system. On a large but finite lattice, however, the behavior of the 
model is similar to the one of a 4-dimensional scalar model: On a finite lattice we regard the 
Kosterlitz-Thouless phase as the "symmetric" phase and the spin- wave phase as the "broken" 



This is explained in the next section 
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phase. In the "broken" phase we have numerically found that there exists a critical value of 
a domain-wall mass rajj which separates a region with a fermionic zero mode on the domain 
wall from one without it. In the "symmetric" phase the critical values of the domain wall 
mass seems to also exist but is very close to its upper bound mo = 1, so that the region with 
a fermionic zero mode is very narrow. Because of the difficulty observed in the numerical 
simulation near m = 1 we cannot exclude a possibility that the existence of the zero mode 
is an artifact of finite lattice size effects. Further simulation we have made on larger lattice 
sizes can not give a definite conclusion. Since, as mentioned before, the phases of the model 
in the infinite volume limit is different from the ones of the 4-dimensional model, we did not 
attempt to increase lattice sizes further, for example 100 2 , to see the fate of zero mode in the 
symmetric phase. Instead we have decided to investigate the (4+l)-dimensional U(l) model 
directly, to obtain the definite conclusion on the existence of zero modes in the symmetric 
phase. 

In this paper, in order to know the fate of the chiral zero mode, we have carried out 
a numerical simulation of a domain-wall model in 4+1 dimensions with a quenched U(l) 
gauge field in the (5 = oo limit. In Sec. 2 , we have defined our domain- wall model with 
dynamical gauge fields. We have calculated a fermion propagator by using a kind of mean- 
field approximation, to show that there is a critical value of the domain-wall mass parameter 
above which the zero mode exist. The value of the critical mass may depend on f3 s , which 
controls the dynamics of the gauge field. In Sec. 3 , we have calculated the fermion spectrum 
numerically using quenched approximation at f3 s = 0.29, 0.5 and at various values of domain- 
wall masses. We have found that in the broken phase (f3 s = 0.5) there exists the range of a 
domain-wall mass parameter in which the chiral zero mode survives on the domain-wall. In 
the symmetric phase (f3 s = 0.29), however, from data on several lattice sizes we have found 
an numerical evidence that the chiral zero mode disappears in the infinite volume limit of 
4-dimensional Euclidean space-time. Our conclusions and some discussions are given in Sec. 
4. 
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II. DOMAIN- WALL MODEL 



A. Definition of the model 

We consider a vector gauge theory in D = 2k + 1 dimensions with a domain- wall mass 
term, which has a kink-like mass term in the coordinate of an extra dimension. This domain- 
wall model is originally proposed by Kaplan ||, and a fermionic part of the action is refor- 
mulated by Narayanan- Neuberger M, in terms of a 2A;-dimensional theory. The model is 
defined by the action 



S — Sn + S, 



(1) 



where Sq is the action of a dynamical gauge field , Sp is the fermionic action. Sq is given 
by 



s g = p E E^-^Mn-^ + AEEO-^IMn.*)]}. (2) 

n,ii>v s n,fi s 

where /x, v run from 1 to 2k , n is a 2fc-dimensional lattice point , and s is a coordinate of 
an extra dimension, [/^(n, s) is a 2fc-dimensional plaquette and C/ M o(n, s) is a plaquette 
containing two link variables in the extra direction. /3 is the inverse gauge coupling for the 
plaquette and (3 S is the one for the plaquette U^d ■ In general , (5 /3 S . The fermion 
action Sf on the Euclidean lattice, in terms of the 2/c-dimensional notation, is given by 

^f = ^EE 4>s{n)i^ [U S)tM (n)ip s (n + n) - U\^{n - n)ip s {n - n) 

n[i s 

t 



+ EEiW MqP r + MZP l Mn 



n s,t 



+ nEE^(") \UsA n )^sin + fi) + U\^{n - n)ip s (n - y) - 2ip s (n) 



(3) 



n/i s 



where s, t are an extra coordinates , Pr/l = |(1 ± 72fc+i) , 

(M ) s>t = U S)D (n)5 s+1)t - a(s)5 S)t 
(Mi) 8tt = U\_ x D (n)6 a -i it - a(s)S s>t . 



(4) 
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Here U s ^(n), U S} d(ti) (D = 2k + 1) are link variables connecting a site (n, s) to {n + ji, s) or 
(n, s + 1), respectively. Because of a periodic boundary condition in the extra dimension , 
s, t run from — L s to L s — 1 , and a(s) is given by 



a(sj = 1 — m sign 



(« + -) sign(L s - s - -) 



1 - too (-5 < s < L s - |) 
1 + m (-L s - ~ < s < -~) , 



(5) 



where mo is the height of the domain-wall mass. It is easy to check that the above fermionic 
action is identical to the one in 2 k + 1 dimensions, proposed by Kaplan HQ . 

In weak coupling limit of both j3 and f3 s , it has been shown that at < m < 1 a desired 
chiral zero mode appears on a domain wall (s = plane) without unwanted doublers. 
Because of the periodic boundary condition in the extra dimension, however, a zero mode 
of the opposite chirality to the one on the domain wall appears on the ant i- domain- wall 
(s = L s — 1). Overlap between two zero modes decreases exponentially at large L s . A 
free fermion propagator is easily calculated and an effective action of a (2+1)- and (4+1)- 
dimensional model including the gauge anomaly and the Chern-Simons term can be obtained 
for smooth background gauge fields J3|||. 

The original Kaplan's do main- wall models in the 4+1 dimension, however, have not 
been investigated yet non-perturbatively, except f3 s = PU| , P2] , P5| and (2+l)-dimensional 
case 



24j. Main question is whether the chiral zero mode survives in the presence of rough 
gauge fields mentioned in the introduction. To answer this question we will analyze the fate 
of the chiral zero mode in the weak coupling limit for j3. In this limit, the gauge field action 
So is reduced to 

S^ftEEl 1 " ReTr [V(n, s)yt(„ + ^ s )] } , (6) 

s n,/i 

where the link variable U Sj o(n) in the extra direction is regarded as a site variable V(n, s)(= 
U St £)(n)). This action is identical to the one of a 2fc-dimensional spin model and s is regarded 
as an independent flavor. The action eq.(^j) is invariant under 



V(n,s)^g(s)V(n,s)g\s + l) , (g(s) 6 G), (7) 

where G is the gauge group of the original model. Therefore the total symmetry of the 
model is G 2Ls , where 2L S , the size of the extra dimension, is regarded as the number of 
independent flavors. We use this (reduced) model for our numerical investigation. 



B. Mean field approximation for fermion propagators 

When the dynamical gauge fields are added even on the extra dimension only, it is 
difficult to calculate the fermion propagator analytically. Instead of calculating the fermion 
propagator exactly, we use a mean-field approximation to see an effect of the dynamical 
gauge field qualitatively. The mean-field approximation we adopt is that the link variables 
are replaced as 



V(n,s) (=U s>D (n)) — > z, 



(8) 



where z is a (n, s)-independent constant. From eq.(^|) the fermion action in a 2/c-dimensional 
momentum space becomes 



S F ^J2M~P) (E^sin(^)5 Sit + [M(z)P R + M\z)P l ] s \ Mp) 



(9) 



s,t,p 



D-l 



{M(z)) s>t = (M (z)) s>t + ^-5 s>t , (M\z)) s , t = {M^z)) s , t + ^-6 s , t , V(p) = £ 2(cosp M - 1), 



(M (z)) s>t = z5 s+ljt - a(s)6 s>t , (M^(z)) s>t = zS B - l t - a(s)8 s>t 



(10) 



(11) 



Following Refs. H,|5]|| especially Ref. 0, it is easy to obtain a mean field fermion propagator 
on a finite lattice with the periodic boundary condition: 



G(p) t 



iJ2^ + M(z)^G L (p)j p L+ ^(-iJ2^ + M\z)^G R (p)} P R 



;i2) 



G L (p) = 



, G R (p) 



(13) 



p 2 + M^(z)M(z) ' p 2 + M(z)M\z) ' 

with p^ = sin(p M ). For large L s where we neglect terms of 0(e~ cLa ) with c > 0, Gl and G# 
are given by 



[Gl(p)] 



s,t 



Be -a+\s-t\ + _ e -a+(«+t) + _ e -a+(2L.- S -t) > ( s? £ > Q) 

A L e- Q + s+Q -* + ^ i?e -«+(i s -s)-a-(L s +t) ? ( s > o, t < 0) 

A L e a - s - a +* + A R e- a - < ~ Ls+s ^- a +^- t \ (s < 0, t > 0) 

C e -a_| a -t| + _ C ) e a_( S +t) + ( Ar _ Q e - a _(2L s +s+t) ? ^ £ < Q) 



(14) 



[Gr{p)U = 



Be -a+\s-t\ + _ e - a+ ( S +t+2) + _ e -a+(2L s -s-t-2) ? ^ t > _l) 

^ e -a+(s+l)+a_(t+l) _|_ ^ e -a+(L s -s-l)-a_(L s +t+l)^ ^ > < _1) 

^ ij , e a -( s + 1 )-«+(*+ 1 ) -|_ ALe -a-(L s +s+l)-a+(L s -t-l) ^ ^ < — 1,£ > — 1) 

Ce -a-\s-t\ + _ C ) e a_( S +t+2) + _ C ) e -a_(2L s+S+t+ 2) ) ^ j < _^ 

(15) 



where 



V(p) 

a± = ^(1 — =f m ) = zfc±, 



a± = arccosh 
A L 
B -- 



p 2 + z 2 + 6 2 



± 



2zb 



± 



1 



1 



a + e 



-Q+ 



c 



1 



(16) 
(17) 
(18) 
(19) 



2a + sinh a + 2a_ sinh a_ 

Like a free fermion theory, the terms of Ar, B and C have no singularity for all z as p — > 
A behavior of A L is, however, different. As p — > A L behaves as 

1 



[(1 - moY ~ A + 0(p 2 ) 

±m 2 n -\(z 2 -l)-ml\ 2 
4m z 2 p 2 



, (0 < m < 1 — z) 



(20) 



, (1 — z < mo < 1). 



A critical value of the domain-wall mass that separates a region with a zero mode and a 
region without zero modes is mjj — 1 — z. Since Al term dominates for 1 — z < m < 1 in 
the Gl (eq. (|I4D) and Gr (eq.(|15D), a right-handed zero mode appears in the s = plane , 
and a left-handed zero mode in the s = L s — 1 plane. For < mo < 1 — z the right- and 
left-handed fermions are massive in all s-planes. Since the terms of Al(Ar) an d B(C) are 
almost same value in this region of mo, a translational invariant term dominates in Gl and 
Gr in the positive (negative) s-layer , so that the spectrum becomes vector-like. 

If z — > 1 , the model becomes a free theory. The propagator obtained in this section 
agrees with the one obtained in Ref. In the opposite limit that z — > , since there is 
no hopping term to the neighboring layers, this model becomes the one analyzed in Ref. 
[ 22[ in the case of the strong coupling limit (3 S = , and in Ref. p3] , in the case that z is 



identified to the vacuum expectation value of the link variables. This consideration suggests 
that the region where the zero modes exist become smaller and smaller as z {1 — z < mo < 1) 
approaches zero. 

What corresponds to z ? Boundary conditions which z satisfies are z = 1 at (3 S = oo 
and z = at (3 S = 0. As explained before the gauge field action of our model is identical to 
that of the U(l) spin system in 4-dimensions. So the most naive candidate El 



is 



z=(V(n,s)), (21) 

which is not invariant under the symmetry (0). In our simulations on L 3 x L 4 x 2L S lattices, 
as an order parameter, we take a vacuum expectation value of link variable calculated with 
rotational technique: 

1 



-E 



5>(n,a) 



) ■ (22) 



L 4 L 3 

Although v defined in eq. (p2"D is always non-zero on finite lattices, it becomes zero in the 
symmetric phase but stays non-zero in the broken phase in the infinite volume limit. Fig. 
|T](a) shows that, even on finite lattices, v behaves as if it was an order parameter: it is 
very small and decreases as the volume increases in the would-be symmetric phase. If the 
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identification that z = v is true, zero modes disappear in the symmetric phase, where v = 0. 
The other choice, which is invariant under (|7|), is the vacuum expectation value of plaquette: 

z 2 = (TrRe{y(n, s)V\n + //, s)}). (23) 

In this identification the zero modes always exist in both phases, since (TrRe{V(n, s)V J< (n + 
//, s)}) is non-zero for all (3 S except (3 S = and is insensitive to which phase we are in, as 
shown in Fig. [1](b). 

III. NUMERICAL STUDY OF (4+l)-DIMENSIONAL U(l) MODEL 
A. Method of numerical calculations 

In this section we numerically study the domain- wall model in 4+1 dimensions with a 
U(l) dynamical gauge field in the extra dimension. As seen from eq.(|^) , the gauge field 
action can be identified with a 4-dimensional U(l) spin model (with 2L S copies). 

Our numerical simulation has been carried out by the quenched approximation. Config- 
urations of U(l) dynamical gauge field are generated and fermion propagators are calculated 
on these configurations. The obtained fermion propagators are gauge non-invariant in gen- 
eral under the symmetry (0). The fermion propagator G(p) Stt becomes "invariant" if and 
only if s = t. Thus, we take the s — s layer as propagating plane (~ "physical space"), and 
investigate the behavior of the fermion propagator in this layer. 

To study the fermion spectrum, we first extract Gl and Gr, defined in eq.(|TI]), from 
the fermion propagator. Assuming eq. (|T3|) , we then obtain corresponding "fermion masses" 
from G~[ 1 (p) and G] i 1 (p) by fitting them linearly in p 2 as follows. 

Gj} =p 2 + M^M ^m) (Right) , (p -> 0), (24) 
Gx 1 =p 2 + MM^ -*m) (Left) , (p -> 0). (25) 

We take the following setup for 4-dimensional momenta. A periodic boundary condition 
is taken for the 1st- 2nd- 3rd-direction and the momenta in these directions are fixed on 
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PiiP2,P3 — 0. An anti-periodic boundary condition is taken for the 4th-direction and the 
momentum in this direction is variable, p 4 = (2n + l)7r/L 4 , n = — L 4 /2, L 4 /2 — 1, where 
L 4 is the number of site of 4th direction. (Our numerical simulations have been performed 
always for L 4 = 32.) 

B. Simulation parameters 

From Fig. |l|(a) it is inferred that the system is in the symmetric phase at (3 S = 0.29 and 
in the broken phase at f3 s = 0.5. Our simulation is performed in the quenched approximation 
on L 3 x 32 x 2L S lattices with L =4, 6, 8 and L s = 8 at (3 S = 0.29 (symmetric phase) , 
where L is a lattice size of 1st- 2nd- 3rd-direction, and on L 3 x 32 x 2L S lattices with L =4, 
6, 8 and L s = 8 at (3 S = 0.5 (broken phase). The coordinate s in the extra dimension runs 
—8 < s < 7. Gauge configurations are generated by the 5-hit Metropolis algorithm. For the 
thermalization first 5000 sweeps are discarded. 

The fermion propagators are calculated by the conjugate gradient method on 50 config- 
urations separated by 100 sweeps. We take the domain-wall mass m = 0.7 , 0.8 , 0.9 , 0.95 
, 0.99 at f3 s = 0.29 and m = 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.6 , 0.9 at f3 s = 0.5. As mentioned 
before, the boundary conditions in 1st- 2nd- 3rd- and 5th-directions are periodic and the 
one in 4th-direction is anti-periodic. Wilson parameter r has been set to r = 1. The fermion 
propagators have been investigated mainly at s — , — 1. These s are the layers where we 
put sources. The layer at s = is the domain wall. Errors are all estimated by the jack-knife 
method with unit bin size. 

C. Fermion spectrum in the broken phase 

The system is in broken phase at f3 s = 0.5. We first consider the fermion spectrum on 
the layer at s — 0. Let us show Fig. |2|, which is a plot of the G~[} and as a function of 
p\ = sin 2 (p 4 ) at m =0.1 and 0.6. (Note we always set Pi,P2,P3 = 0.) In the limit p 4 — > 0, 
Gft 1 remains non-zero at both mo, while G~j} vanishes at mo = 0.6. We obtain the value 
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of m 2 f, which can be regarded as the mass square in 4-dimensional world, by the linear fit 
in p\ near p\ = 0, and plot mj as a function of mo in Fig. [| The mass of right-handed 
fermion, obtained from G^ 1 , becomes very small (less than 0.1) at mo larger than 0.35, so 
we conclude that the critical value is ~ 0.35. Whenever the domain-wall mass is larger 
than this value , this model produces the right-handed chiral zero mode on the domain 
wall at s = 0. From these results above we conclude that the domain-wall model with the 
dynamical gauge field on the extra dimension (i.e. the weak coupling limit of the original 
model) can survive the chiral zero mode on the domain-wall, at least deep in the broken 
phase. 

Also in Fig. ^|, solid lines stand for the inverse propagators obtained from the mean-field 
propagator with appropriately tuned parameter z, and the lines show that the behavior of 
the fermion propagator is well described by the mean- field propagator. 

An important question here is what a fermion spectrum is in the scaling limit. If the 
fermion spectrum stays chiral in the limit, it should stay chiral also in the symmetric phase, 
since the phase transition is continuous. This means that, in order to determine the fermion 
spectrum in the scaling limit even from the broken phase, we have to know the spectrum in 
the symmetric phase. Therefore, from the knowledge of the fermion spectrum obtained in 
the broken phase so far, we cannot draw any conclusions on the fermion spectrum, chiral or 
vector-like, in the scaling limit. 

D. Fermion spectrum in symmetric phase 

The system is in the symmetric phase at (3 S = 0.29. The fermion propagator is analyzed 
in the same way as in the broken phase. 

In Fig. |], we have plotted mass m./ of the right- and left-handed modes at s = as a 
function of itlq. On a 4 3 x 32 x 16 lattice, it seems that the right-handed fermion becomes 
massless at mo larger than 0.95, while the left-handed fermion stays massive at all mo, so 
that the fermion spectrum on the domain wall is chiral. However as the lattice sizes become 
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larger, for example 6 3 x 32 x 16 and 8 3 x 32 x 16, mass differences between the right- and 
the left-handed modes become smaller. This suggests that the fermion spectrum becomes 
vector-like in the infinite volume limit. From this data alone, however, we cannot exclude 
a possibility that the critical mass m§ is very close to 1.0, since the fermion mass near 
m = 1.0 is very small. In order to make a definite conclusion on the absence of chiral zero 
modes in the symmetric phase, we try to fit G]~} and G^, 1 at given mo using the form of 
mean-field propagator, eqs.fljl]) and fllED , with the fitting parameter z. We show the quality 
of this fit in Figs. ||] (a) and (b). These figures show that the fermion propagator is well 
described by the mean-field propagator if the fitting parameter z is chosen such that x 2 is 
minimized. We then have plotted z obtained by the fit as a function of 1/L in Figs. |6| (a) 
and (b), where L's take 4, 6, 8 on the L 3 x 32 x 16 lattices. The values of z's are almost 
independent of m at the each 1/L except the right-handed modes at m = 0.99. The solid 



circles represent the order parameter v defined in eq. (|22|) . The behaviors of z's at different 
m are almost identical each other and are very similar to that of v except the right-handed 
ones at m = 0.99. We think that the observed deviation of z (right-handed) at m = 0.99 
from those at other m 's is not a real effect but a statistical fluctuation, since the \ 2 for the 
fit at m = 0.99 is almost flat in the region between z = and z ~ 0.3. Ignoring this data 
at mo = 0.99 we see that z can be identified with v, and therefore z becomes zero as the 
lattice size goes to infinity since the order parameter v defined in eq . (|22|) becomes zero in 



the infinite volume limit. This analysis leads us to a conclusion that the fermion spectrum 
in the symmetric phase is vector-like. 

Fig. |] also shows that the fermion masses tend to approach the line rrif = 1 — m , which 
is the value of the mass for free Wilson fermion at a positive s-layer, as the lattice volume 
increases. Furthermore in Fig. [7] as a function of mo, we have also plotted rrif at s — — 1 
(a negative s-layer). As the volume increases the mass difference between right-handed 
and left-handed fermions becomes smaller and both masses seem to approach to the line 
rrif = 1 + m , the value of the mass for free Wilson fermion at a negative s-layer. This data 
also supports our conclusion on the absence of chiral zero modes in the symmetric phase. 
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In summary our results of m/ at both positive s- and negative s-layer suggest that 
chiral zero modes in the symmetric phase disappear in the infinite volume limit. This 
concludes that the original Kaplan's model fails to describe lattice chiral gauge theories in 
the symmetric phase. 



IV. CONCLUSIONS AND DISCUSSIONS 

Using the quenched approximation, we have performed the numerical study of the 
domain-wall model in 4+1 dimensions with the U(l) dynamical gauge field on the extra 
dimension. From this study we obtain the following results. In the broken phase of the 
gauge field, there exists the critical value of the domain-wall mass separating the region 
with a chiral zero mode and the region without it in the same case of (2+l)-dimensional 
model. At the domain-wall mass larger than its critical value, a zero mode with one chi- 
rality exists on the domain wall. In the symmetric phase, on the other hand, our data on 
L 3 x 32 x 16 with L =4 ,6,8 lattices suggest the absence of chiral zero modes in the infinite 
volume limit, though the chiral zero mode seems to exist on finite lattices. We also found 
that fermion propagators obtained through numerical simulations on finite lattices are well 
described by the mean-field propagator with z ~ v. Since the existence of chiral zero modes 
in the symmetric phase is essential for the success of the original domain-wall model, our 
results for the (4+l)-dimensional model indicate that the original domain- wall model cannot 
work as lattice chiral gauge theories. 

In 2+1 dimensions we could not conclude whether the fermionic zero mode exists or 
not because of the difficulty for the simulation near m = 1 and the absence of the order 
parameter in the infinite volume limit for the 2-dimensional U(l) model. We try here to 
extract the parameter z of the (2+l)-dimensional model, by the method used for the (4+1)- 
dimensional model. In Figs, p] (a) and (b), the parameter z is plotted as a function of 1/L. 



The behavior of z is almost identical to that of order parameter v in eq.(22) but not to the 



square root of plaquette in eq. (|2"B"D . This new analysis shows that zero modes are absent in 
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the symmetric phase in the (2+l)-dimensional U(l) model as well as the (4+l)-dimensional 
model. Therefore we conclude that U(l) chiral gauge theories can not be constructed via 
an original domain-wall model, regardless of its dimensionality. 

One of the remaining question is whether the above negative conclusion also holds for 
other gauge groups such as SU(2). The answer is not so straightforward: For example, the 
2-dimensional SU(2) spin model has a symmetric phase only. We wonder whether chiral zero 
modes are absent in the symmetric phase even if the gauge field of the (2+1) -dimensional 
domain- wall model becomes "smooth" for large but finite f3 s . To answer this question we 
are currently investigating the (2+l)-dimensional SU(2) domain- wall model. 
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FIGURES 
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FIG. 1. (a) Vacuum expectation value of link variables v and (b) Vacuum expectation value 
of plaquette w on L 3 x 32 x 16 lattices with L =4(circles) , 6(diamonds) , 8(triangles) as a function 
of/3,. 
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Fig. 2(a) 

ps=0.5 , 4 3 x32x16 , s=0 
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Fig. 2(b) 
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FIG. 2. (a) [GrJoo as a function of sin 2 (j>4) and (b) [GlJoo as a function of sin 2 (j>4) with 
PiiP2,P3 = at (3 S = 0.5 on a 4 3 x 32 x 16 lattice, at mo=0.1 (solid circles) and 0.6 (solid 
triangles). Solid lines of both figures stand for the ones obtained from the mean- field propagator 
with the fitted parameter z. 
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(3s=0.5 , L 3 x32x16 , s=0 
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FIG. 3. rrif vs. mo at (3 S = 0.5 on a L 3 x 32 x 16 lattices with L =4(circles) , 6 (up triangles) 
and 8 (down triangles) in the case of putting a source on the domain wall at s = 0, for the 
right-handed fermion (solid symbols) and the left-handed fermion (open symbols). 
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Fig. 4 

Ps=0.29 , L 3 x32x16 , s=0 
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FIG. 4. nif vs. m at j3 s = 0.29 on L 3 x 32 x 16 lattices with L =4 (circles) , 6 (up 
triangles) and 8 (down triangles) in the case of putting a source on the domain wall at s = 0, for 
the right-handed fermion (solid symbols) and the left-handed fermion (open symbols). Solid line 
corresponds to 1 — mo. 
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Fig. 5(b) 

ps=0.29 , 4 3 x32x16 , S=0 
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FIG. 5. (a) [G R ] oJ as a function of sin 2 (j?4) and (b) [Gl] q as a function of sin 2 (p4) with 
Pi,P2,P3 = at (3 S = 0.29 on a 4 3 x 32 x 16 lattice, at mo=0.7. Solid lines of both figures stand 
for the ones obtained from the mean- field propagator with the fitted parameter z. 
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Fig. 6(b) 

Ps=0.29 , L 3 x32x16 , S=0 
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FIG. 6. (a) z (right-handed) and i; vs. 1/L at (3 S = 0.29 and (b) z (left-handed) and v vs. 
1/L at 0s = 0.29. on L 3 x 32 x 16 lattices with L =4, 6 and 8 at m = 0.7 (up triangles), 0.95 
(diamonds) and 0.99 (down triangles), in the case of putting a source on the domain wall at s = 0. 
Solid circles stand for the vacuum expectation value of link variable (:the order parameter). 



24 



m, 



Fig. 7 

ps=0.29 , L 3 x32x16 , 50 conf. , s=-1 
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FIG. 7. mf vs. m at [3 S = 0.29 on L 3 x 32 x 16 lattices with L =4 (circles) , 6 (up triangles) 
and 8 (down triangles) in the case of putting a source at s = —1, for the right-handed fermion 
(solid symbols) and the left-handed fermion (open symbols). Solid line corresponds to 1 + m®. 
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Fig. 8(a) 

L 2 x32 , ps=0.5 (in the symmetric phase) , s=0 
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Fig. 8(b) 

L 2 x32 , ps=0.5 (in the symmetric phase) , s=0 
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FIG. 8. (a) z(right-handed) and v vs. 1/L at (3 S = 0.5 and (b) z (left-handed) and v vs. 1/L 
at p s = 0.5 in the symmetric phase on I? x 32 lattices with L =24 and 32 at mo = 0.7 (squares), 
0.8 (up triangles), 0.9 (diamonds) and 0.99 (down triangles), in the case of putting a source on the 
domain wall at s = 0. Solid circles stand for the vacuum expectation value of link variable (:the 
order parameter). 
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